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Abstract 

Rationale: The degree to which tuberculosis (TB) is transmitted between persons is variable. Identifying the factors that 
contribute to transmission could provide new opportunities for TB control. Transmission is influenced by host, bacterial and 
environmental factors. However, distinguishing their individual effects is problematic because measures of disease severity 
are tightly correlated, and assessing the virulence of Mycobacterium tuberculosis isolates is complicated by epidemiological 
and clinical confounders. 

Objectives: To overcome these problems, we investigated factors potentially associated with TB transmission within 
households. 

Methods: We evaluated patients with smear-positive (S2+), pulmonary TB and classified M. tuberculosis strains into single 
nucleotide polymorphism genetic cluster groups (SCG). We recorded index case, household contact, and environmental 
characteristics and tested contacts with tuberculin skin test (TST) and interferon-gamma release assay. Households were 
classified as high (>70% of contacts with TST>10 mm) and low (<40%) transmission. We used logistic regression to 
determine independent predictors. 

Result:Fmm March 2008 to June 2012, we screened 293 TB patients to enroll 124 index cases and their 731 contacts. There 
were 23 low and 73 high transmission households. Index case factors associated with high transmission were severity of 
cough as measured by a visual analog cough scale (VACS) and the Leicester Cough Questionnaire (LCQ), and cavitation on 
chest radiograph. SCG 3b strains tended to be more prevalent in low (27.3%) than in high (12.5%) transmission households 
(p = 0.1 1). In adjusted models, only VACS (p<0.001) remained significant. SCG was associated with bilateral disease on chest 
radiograph (p = 0.002) and marginally associated with LCQ sores (p = 0.058), with group 3b patients having weaker cough. 

Conclusions: We found differential transmission among otherwise clinically similar patients with advanced TB disease. We 
propose that distinct strains may cause differing patterns of cough strength and cavitation in the host leading to diverging 
infectiousness. Larger studies are needed to verify this hypothesis. 
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Introduction 

Tuberculosis (TB) remains one of the most important threats to 
global health. Although TB is the prototypic disease with airborne 
transmission [1], there is both experimental [2,3] and epidemio- 
logic [4,5,6,7,8,9] evidence of marked variability in transmission of 
M. tubercubsis (MTB) from patients with pulmonary TB. The 
reasons underlying this variability remain poorly understood in 



part because of the inherent complexity of studying a biological 
phenomenon resulting from a complex web of interactions 
between two human hosts-the source case and the exposed 
contact-and the infecting pathogen within an array of environ- 
ments. Yet, these differences have important implications for 
vaccine development and global TB control. 



PLOS ONE | www.plosone.org 



1 



July 2014 | Volume 9 | Issue 7 | e1 00984 



M. tuberculosis Transmission within Households 



Most studies of transmission and disease in household contacts 
have focused on host and environmental factors [7,10]. These 
studies fail to explain differences in epidemiologic behavior of 
clinical isolates of MTB reported in outbreaks. In fact, there is a 
growing body of experimental evidence suggesting that MTB 
strains have variable degrees of virulence 
[11,12,13,14,15,16,17,18,19]. Although a great deal is known 
about the genome and proteome of MTB, the genetic basis for 
most of the cited differences in strain behavior is unknown in most 
cases [20]. Despite this growing body of literature, there have been 
no systematic studies designed to investigate the role of strain 
difference as a determinant of the critical epidemiologic and 
clinical parameters such as transmissibility, resistance to infection 
and progression to disease. Importantly, most of the cited studies 
investigating the virulence potential of MTB isolates have used 
banked clinical isolates with limited knowledge about the 
epidemiologic and clinical conditions present at the time the 
isolates were collected [11,12,13,14,15,16,17,18,19,20]. 

We conducted a household contact study in Vitoria, Brazil to 
test the hypothesis that genetic strain variation in MTB was a 
determinant of the extent of transmission of infection in household 
contacts of smear-positive cases of pulmonary TB. We sought to 
determine whether some clinical strains of MTB were more 
transmissible than others, while controlling for other host and 
environmental factors involved in transmission. 

Methods 

Ethics Statement 

The study was approved by the Institutional Review Boards of 
Boston University Medical Campus and New Jersey Medical 
School - Rutgers University (formerly University of Medicine and 
Dentistry of New Jersey), the Comite de Etica em Pesquisa do 
Centro de Ciencias da Saude - Universidade Federal do Espirito 
Santo and the Comissao Nacional de Etica em Pesquisa (CONEP). 
We obtained written informed consent and assent in Portuguese in 
accordance with age-specific ethical guidelines of participating 
institutions. 

Study population 

This study was conducted at the The Nucleo de Doencas 
Infecciosas (NDI) located in Vitoria, the capital city of the State of 
Espirito Santo, Brazil. The NDI has organized a network of five 
laboratories in the metropolitan region of Vitoria that serve 16 TB 
clinics distributed in the municipalities of Vitoria, Cariacica, Serra 
and Vila Velha. With approximately 1 ,400 cases/ year, the annual 
TB incidence in Espirito Santo is 38/100,000 inhabitants. The 
prevalence of HIV infection in Espirito Santo is < 1 % in the 
general population, and 7% in TB cases [21]. 

All consecutive pulmonary TB patients identified through the 
NDI clinic network were eligible to participate in this study, 
provided they fulfilled the following inclusion criteria: 1) age £18 
years with cough >3 weeks; 2) new TB episode with >1 sputum 
specimen with acid-fast bacilli (AFB) >2+ with subsequent MTB 
growth in culture; and 3) living with S3 household contacts. We 
excluded index cases who were HIV-infected (or refused HIV 
testing), had a history of TB treatment, or who were too ill to 
consent, unable to understand, or to comply with the study 
protocol. To minimize differences in exposure time between study 
households, index TB cases were screened and enrolled within the 
first 2 weeks after they first presented to the municipal TB clinic. A 
household contact was defined as an individual of any age fulfilling 
at least one of the following criteria of close contact with the index 
TB cases for >3 months before enrollment: 1) sleeping under the 



same roof S5 days/week; 2) sharing meals £5 days/week; 3) 
watching TV nights or weekends, and; 4) other significant contact 
(85% of these visited the household >18 days/month). 

Measurements 

TB cases. We collected demographic and clinical informa- 
tion about TB patients and recorded their cough severity using two 
measurements: i) a self-reported visual analog cough scale (VACS) 
[22], and ii) the Leicester Cough Questionnaire (LCQ) [23,24]. 
The LCQ is a 19-item, patient-derived questionnaire that 
measures the physical, psychological and social effects of chronic 
cough in adults; the final score (range 3-21) is inversely related to 
cough severity (e.g. higher scores indicate weaker cough). Because 
of the inherent complexity of the LCQ, study staff administered 
the LCQ only to TB patients deemed reliable historians (over 
95%). We obtained up to three sputa specimens for AFB smear 
microscopy (auramine O fluorescent stain) [25] and culture 
(Ogawa-Kudoh method). MTB isolates underwent single nucleo- 
tide polymorphism (SNP) analysis as described previously [26] and 
were categorized into one of nine SNP genetic cluster groups and 
subgroups (SCG), as described previously [27]. The radiological 
extent of lung disease was graded on a four-category ordinal scale 
(normal, minimal, moderate and far-advanced) by an experienced 
physician [28]. All patients were offered standard TB treatment 
according to Brazilian guidelines [29,30]. Study staff visited the 
participants' dwellings to verify the identity of each contact, 
measure individual contact time, and to perform an environmental 
evaluation (crowding and ventilation). 

Household contacts. We followed recommendations of the 
Brazilian TB Control Program for household contact investiga- 
tions [29,30]. We recorded demographic and clinical character- 
istics of contacts, and evaluated them for MTB infection with both 
TST and IGRA (Quantiferon Gold-In-Tube, Quiagen). To ensure 
accurate TST readings, only staff trained by the National TB 
Program were responsible for TST application, and we completed 
inter- and intra-reader evaluations (kappa >90%) prior to opening 
study enrollment. Briefly, we placed two units of R23 (SSI, 
Denmark) provided by the Brazilian Ministry of Health on the 
forearm of contacts using the Mantoux method; the diameter of 
induration was measured in millimeters between 72 and 96 hours. 
Household contacts with a TST<10 mm at baseline were re- 
tested after 8-12 weeks to identify TST conversion, using two 
different criteria for conversion: Criterion 1 (Brazilian Guidelines) : I st 
TST<10 mm, 2 nd TST>10 mm, and difference between 1 st and 
2 nd TST>10 mm. Criterion 2 : 1 st TST<5 mm, 2 nd TST>10 mm, 
and difference between l ht and 2 nd TST&6 mm [31]. At 8-12 
weeks, blood for IGRA testing was obtained before TST 
placement to reduce TST-induced IGRA boosting [32]. Following 
Brazilian recommendations [29], we referred contacts with TST> 
10 mm and secondary TB suspects for evaluation and treatment 
by the corresponding Municipal TB Clinic. 

Household transmission categories 

We categorized households into MTB transmission categories 
according to the percentage of contacts in each household that had 
a TST>10 mm by 8-12 weeks: low (LT) transmission (£40%) or 
high (HT) transmission (>70%). This definition includes contacts 
with TST>10 mm at either baseline or week 8-12, providing a 
more complete measure of the total infection burden for each 
household. Originally, 25% and 75% cutoffs for LT and HT were 
planned but relaxed to 40% and 70%, respectively, based solely on 
increasing the numbers in the two categories; no other factors were 
examined when the cutoffs were revised. A group classified as 
intermediate transmission (IT) households (41-69% contacts with 
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TSTslO mm) were also enrolled to serve as controls. Contacts 
with a history of previous TB disease (N = 13) and with missing 
results on 1 st or 2 nd TST (N = 44) were retained in the general 
analysis but did not contribute to the determination of the 
transmission category; their inclusion or exclusion had little effect 
on the category (data not shown). Two contacts were diagnosed 
with co-prevalent TB disease and based on symptom history were 
determined to have been infected by the index case and not the 
reverse; these two contacts were considered TST-positive for 
purposes of the determination of the transmission category. 

Statistical methods 

Because this study used a household contact design, character- 
istics are at both a household level (including data about the index 
case) and at an individual contact level, nested within households. 
Households were classified into three phenotypes, i.e., low, 
intermediate, and high transmission. Comparisons were made 
between LT and HT households, and also across the three 
phenotypes. Unadjusted associations of household level data to 
transmission phenotype assumed independence across households. 
Logistic regression was used to obtain odds ratios of household 
level characteristics, but in the case of sparse data, exact logistic 
regression and corresponding odds ratios and tests were substitut- 
ed [33]. Characteristics from individual contacts within a 
household were potentially correlated and, therefore, were 
evaluated using a generalized estimating equation (GEE) approach 
with an independent working correlation structure. We assumed a 
binomial distribution with a logit link when comparing LT and 
HT households and a multinomial distribution with a cumulative 
logit link when comparing LT, IT, and HT households. Multiple 
regression modeling focused on comparing LT vs. HT households 
and used the GEE approach to account for within household 
correlation. An adhoc approach was used to determine charac- 
teristics for the final multiple regression model. An automatic 
stepwise selection procedure was used to select household level 
characteristics, and then individual contact level covariates were 
candidate for addition to the model in a manual stepwise fashion. 
SCG and contact age were then added to the model. Interactions 
of SCG with measures of disease severity used simplified models 
because of estimation problems. Covariates of questionable 
reliability were excluded from consideration. Because of their 
strong negative correlation (Pearson correlation —0.72*) VACS 
and LCQ were not candidates for inclusion in a single 
multivariable model. Direct comparison of VACS and LCQ_ 
showed that VACS was a superior predictor of transmission 
category based on both Akaike Information criterion (AIC) and 
likelihood score statistic. Moreover, LCQ,is time consuming and 
difficult to administer and score. Therefore, multivariable analyses 
only considered inclusion of VACS. All analyses were conducted 
in SAS 9.2 and testing is two-sided. 

Results 

Between February 2008 and June 2012, we identified 1,071 
potentially eligible TB patients, 743 were determined ineligible 
and not referred to the study, 35 refused study participation, and 
293 underwent screening; of these, 230 families were confirmed 
eligible, but 106 were not enrolled (Figure 1). This analysis 
includes 124 index cases and their 731 household contacts: 23 LT 
families (132 household contacts), 73 HT families (416 contacts), 
and 28 IT families (183 contacts) that were selected for enrollment 
as control households. 



Index TB patients and dwellings 

The median age of TB cases was 36 years (range 25, 46), 64% 
were male, and most had advanced TB disease as measured by 
duration of cough (median 12 weeks), sputum AFB smear (82% 
AFB 3+), or extent of lung disease on chest radiograph (Table 1). 
The severity of cough (as measured by VACS and LCQ) and the 
presence of cavitation on chest radiograph in index TB cases 
increased as the proportion of household members with MTB 
infection increased. There were no differences in socio-economic 
variables (marital status, highest education level, occupation and 
average salary) between household transmission categories, except 
that families living in HT households were more likely to rent 
(P = 0.012). When compared to LT households, contacts from IT 
and HT families were more likely to self-report close contact with 
another person with TB other than the index case (P = 0.02). 
Study dwellings were similar in terms of building materials, 
sanitation, size, crowding, ventilation, and other examined factors 
(Table SI). 

Genotypic strain analysis 

The two most prevalent strains among TB patients were group 
5 (73.9%) and group 3b (16.8%) and their distribution across the 
three household transmission categories differed qualitatively, 
although not statistically significantly (omnibus P = 0.38) (Table 1). 
Group 3b strains were more prevalent in LT (27.3%) and IT 
(20.0%) households compared to HT (12.5%) households; 
conversely, group 5 strains were most prevalent in HT (76.4%) 
households compared to IT (68.2%) and LT (72.0%) households. 
The median [interquartile range (IQR)] LCQ score was 12 (7, 21) 
for patients infected with SCG 3b strains, 11 (6, 21) for SCG 5, 
and 10 (5, 16) for other SCG groups combined (p = 0.056), 
trending toward less coughing in patients infected with a SCG 3b 
isolate. Similarly, the median (IQR) VACS was 6 (0, 10) for SCG 
3b TB patients (N = 19), 7 (0, 10) for SCG 5 (N = 83), and 9 (3, 10) 
for all other groups (N=ll) combined (p = 0. 11), indicating 
patients infected with a SCG 3b strain had a trend toward a 
weaker cough (Figure 2). Patients with SCG 3b isolates were less 
likely to have bilateral disease (15.0%), vs. SCG 5 (51.2%) and 
other SCG (72.7%) (p = 0.002), and less likely to have infiltrates in 
the upper lung: SCG 3b (80.0%), 5 (95.3%), and 4 or 6a (100.0%), 
with the differences being marginally significant (p = 0.060). 
Similarly, cavitations were present in SCG 3b (65.0%), 5 
(71.8%), and 4 or 6a (90.0%) but these differences were not 
statistically significant (p = 0.39). 

Household contacts 

This analysis includes a total of 731 household contacts that 
fulfilled the study definition of close and prolonged contact time 
with the index case (categorized by the first definition met): 63% 
slept under the same roof S5 days/week, 12% watched TV on 
nights/weekends, 1 1% shared meals &5 days/week and 14% had 
other significant contact (85% visited the household >18 days/ 
month). Most (90%) households were evaluated within ±14 days 
of the index TB cases starting anti-tuberculous treatment. 
Household contacts were young (median age 20 years), but 
covering a broad age range (3 months to 87 years); 56.3% were 
female; 77.0% had a BCG vaccination scar; 99% had no prior 
diagnosis of MTB and 98% reported no HIV-infection (Table 1). 
Enrolled contacts spent considerable time in close contact with 
their family member with TB. Contacts were comparable across 
transmission categories in terms of socio-economic condition 
(marital status, education level, occupation and salary), family 
composition (i.e. no new contacts left the household or died in 3 
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Potentially Eligible: 



TB suspect to clinic 
N=1071 



Examined for Eligibility: 



Confirmed Eligible: 



Screened 
N=293 



Eligible 
N=230 



Included in Study: 



TB cases 
Enrolled 
N=124 



Analyzed: 



Low 
Transmission 

IC: N=23 
HHC: N=132 



• Not eligible (743) 
•Refused (35) 



• Not eligible (60) 

• Refused (3) 



• Intermediate transmission (63) 

• Negative contaminated culture (21) 

• <3 HHC eligible (16) 

• IC Lives for <3 months with HHC (4) 

• HTV positive (1) 

• Unable to contact (1) 



Household 
contacts 
screened 
N=759 



\ 






Eligible 
N=735 


\ 




f 


Enrolled 
N=731 



- Refused (15) 

• Unable to contact (5) 

• Previous TST positive (3) 

• Previous TB disease (1) 



• Died before enrollment (2) 

• Unable to contact (1) 

• HHC definition notmet (1) 



Intermediate 
Transmission 

IC: N=28 
HHC: N=183 



High 
Transmission 

IC: N=73 
HHC: N=416 



Figure 1. Study profile. 
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months prior to study inception), alcohol and tobacco consump- 
tion, and history of co-morbidities. 

At baseline, TST was positive in 57.3% of contacts and the 
distribution varied according to transmission category (by design) 
and contact age (Table 1). Whereas the percentage of contacts 
with TST>10 mm living in LT and IT households increased with 
age, it was similar across age strata in HT households. Overall, the 
frequency of IGRA positivity at 8-12 weeks was approximately 
15% lower than TST but showed a similar age-dependent increase 
in prevalence of infection. The median TST diameter in contacts 
with TST>10mm (1 st or 2 nd TST) increased in step with 
household transmission categories (P<0.001). 



Household transmission categories 

Transmission categories of households changed considerably 
between baseline and week 8-12 as a result of TST conversion 
(Tables 2 and 3). At baseline, 41, 30 and 53 of the sampled 
households were LT, IT and HT, respectively; by 8-12 weeks, 
43% of LT and IT household progressed and only 23 remained as 
LT households (Table 2). The use of two different criteria to define 
TST conversion identified the same number of contacts with 
newly acquired infection, however the same contacts were not 
identified by both criteria (Table 3); of the 56 contacts that 
underwent TST conversion by either criterion, 44 (79%) fulfilled 
both TST conversion criteria. Most TST conversions occurred in 
households initially classified as low or intermediate households, as 
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Table 1. Characteristics of index tuberculosis cases, household contacts 
tuberculosis transmission category. 


and study dwellings at baseline by Mycobacterium 




Characteristic 1 


Final Household transmission category 




P value 2 (Multinomial) P value 2 (LT vs. HT) 




Total 


Low (LT) 


Intermediate 


High (HT) 




Index Cases 


N = 124 


N = 23 


N = 28 


N = 73 




Age (years) 






3R no A&\ 


33 [22,46] 


0.34K 0.47W 


Sex 










0.77C 0.58C 


Male 


Rn if,A ^l 


1 a i&n q\ 


1 7 (60 7) 


49 (67.1) 




Female 


44 (35 5) 


y py. 1 ) 


1 1 fiq 31 


24 (32.9) 




Cough visual analog scale 


7 [4,9] 


5 [3,6] 


6 [4,9] 


8 [6,9] 


0.003K 0.002W 


Leicester cough questionnaire 


11 [10,14] 


13 [11,17] 


12 [10,16] 


11 [9,13] 


0.007K 0.007W 


Weeks sick prior to enrollment 


12 [4,20] 


8 [4,24] 


8 [4,18] 


1 2 [8 24] 


0.26K 0.57W 


Housing 










0.01 7T 0.1 IE 


Owned 


1 07 (86.3) 


22 (95.7) 


27 (96.4) 


58 (79.5) 




Rented 


17 (13.7) 


1 (4.3) 


1 (3.6) 


1 5 (20.5) 




Cni i-t-i I rr\ APR cmpar i Kin Kf^cl" f~1 r\p^\ 
-J \J U 1 U 1 1 1 r\r D 3 1 1 Itral 1, 1 II U 1 1 t;:> I UlaUCJ 










o 4Qr o 3nr 


2+ 


22 (17.7) 


6 (26.1) 


4 (14.3) 


12 (16 4) 




3+ 


102 (82.3) 


17 (73.9) 


24 (85.7) 


61 (83.6) 




Extent of lung disease on chest 
radiograph 










0.20J 0.24T 


Normal 


2. {]./) 


1 IA 31 


i rs 7i 

1 \3. /} 


0 (0.0) 




Minimal 


7 fK R\ 

/ \D.O) 


1 IA 31 


i n a\ 


4 (5.6) 




Moderate 


63 (52 1) 


1 3 (56 5) 


15 (55 6) 


35 (49.3) 




Far advanced 


AQ 1AH ^\ 


8 (34 8} 


Q /33 31 


32 (45.1) 




Cavitations 










0.01 2T 0.023T 


Absent 


35 (28.9) 


10 (43.5) 


11 (40.7) 


14 (19.7) 




Present 


86 (71.1) 


13 (56.5) 


16 (59.3) 


57 (80 3) 




^MP t~\\ nrm in 
jnr 1 U :> It 1 y 1 yJ U [J 










0 3RT 0 46F 

U.JO 1 u.wn 


3b 


20 (16.8) 


6 (27.3) 


5 (20.0) 


Q n 2 51 




4 


2 (1.7) 


0 (0.0) 


1 (4.0) 


1 n 4i 




5 


88 (73.9) 


15 (68.2) 


18 (72.0) 


55 (76 4) 




6a 


8 (6.7) 


1 (4.5) 


1 (4.0) 


6 (R 31 




5 or 6a 3 


1 (0.8) 


0 (0.0) 


0 (0.0) 


1 (1 41 




Number of contacts enrolled in 
study 


5 [4,7] 


5 [4,7] 


6 [5,7.5] 


5 [4,7] 


0.1 OK 0.98W 


Household Contacts 


N = 731 


N= 132 


N = 183 


N = 416 




Age (years) 


20 [9,39] 


18.5 [8.5,38.5] 


25 [12,41] 


19 [9,37] 


0.65S 0.80S 












0.85S 0.40S 


Male 


325 (44.5) 


63 (47.7) 


76 (41 .5) 


1 86 (44.7) 




Female 


406 (55.5) 


69 (52.3) 


107 (58.5) 


230 (55.3) 




BCG scar 










0.72S 0.86S 


Present 


563 (77.0) 


107 (81.1) 


130 (71.0) 


326 (78.4) 




Absent 


137 (18.7) 


21 (15.9) 


40 (21.9) 


76 (18.3) 




Uncertain 


31 (4.2) 


4 (3.0) 


13 (7.1) 


14 (3.4) 




Close contact with person with 
TB other than index case 










0.020S <0.001S 


Yes 


159 (21.8) 


7 (5.3) 


46 (25.1) 


106 (25.5) 




No 


570 (78.1) 


124 (93.9) 


137 (74.9) 


309 (74.5) 




TST and IGRA status, n/N (%) 


1 st TSTS10 mm (at baseline) 


407/710 (57.3) 


17/132 (12.9) 


78/180 (43.3) 


312/398 (78.4) 




1 st or 2 nd TSTalO mm 
(at 8-12 weeks) 


488/691 (71.7) 


27/127 (21.3) 


103/174 (59.2) 


358/380 (94.2) 




<5 years 


48/85 (56.5) 


0/20 (0.0) 


9/23 (39.1) 


39/42 (92.9) 
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Characteristic 1 


Final Household transmission category 


P value 2 (Multinomial) P value 2 (LT vs. HT) 




Total 


Low (LT) 


Intermediate 


High (HT) 


6-15 years 


129/190 (67.9) 


3/36 (8.3) 


14/36 (38.9) 


112/118 (94.9) - 


^16 years 


311/406 (76.6) 


24/71 (33.8) 


80/115 (69.6) 


207/220(94.1) - 


Median diameter (mm) 


All 


14 [7,17] 


14 [7,17] 


10 [2,15] 


16 [13,18] 


Among TST>10 mm (only) 


16 [13,19] 


14 [12,17] 


15 [11,19] 


16 [14,19] <0.001S <0.001S 


Positive IGRA (at 8-12 weeks) 


322/566 (56.9) 


18/101 (17.8) 


62/149 (41.6) 


242/316 (76.6) - 


<5 years 


19/47 (40.4) 


0/10 (0.0) 


3/13 (23.1) 


1 6/24 (66.7) 


6-15 years 


95/160 (59.4) 


5/30 (16.7) 


9/30 (30.0) 


81/100(81.0) - 


>16 years 


208/359 (57.9) 


13/61 (21.3) 


50/106 (47.2) 


145/192 (75.5) - 



'n (%) or median [IQR], unless otherwise specified. 

2 C = Chi-squared, E = Exact, K = Kruskal-Wallis, J = Jonckheere-Trepstra, T = Cochran-Armitage Trend, M = McNemar's Test, S = Score test from generalized estimating 
equation (GEE) logistic model; W = Wilcoxon rank sum test. 

3 The SNP cluster group for one index case could not be distinguished between 5 or 6a. This participant's result is shown in the table but is excluded when performing 
statistical tests. 

Missing data: Leicester cough questionnaire (6), cough visual analog scale (6), chest radiograph (3), contact with other TB person than index case (1 missing, 1 don't 
know), 1 st TST (21), 1 st or 2 nd TST (40), IGRA (161 missing, 4 indeterminate), SNP (5). 
Infection status of contacts is shown at baseline and at study conclusion (8 weeks). 
doi:1 0.1 371 /journal.pone.01 00984.t001 



the proportion of contacts "at risk" of TST conversion living in 
HT households was lower by definition. Adolescent contacts were 
more likely to convert their TST than adults (conversion criterion 
1, P = 0.025). 

Factors associated with high MTB transmission 

We analyzed index TB patient, household contact and 
environmental factors associated with HT households (Table 4). 
In unadjusted analyses, factors associated with HT were the 
presence of cavitation in chest radiograph [Odds Ratio (OR) 3.13, 
95% confidence interval (CI): 1.13, 8.69; P = 0.028] and stronger 
cough as measured by VACS (OR 1.36 per unit increase, 95% CI: 
1.13, 1.67; P = 0.001) or LCQ,(OR 1.24 per unit decrease, 95% 
CI: 1.07, 1.45; P= 0.005) in index TB cases. A group 5 SCO MTB 
isolate had more than twice the estimated odds of being associated 
with HT compared to group 3b, although it was not statistically 
significant (OR 2.44, 95% CI: 0.72, 7.93; P = 0.14). In a logistic 
regression model adjusted for other covariates with a GEE 
approach for estimation, the only variable that remained 
independentiy associated with HT was severity of cough in the 
index TB case (VACS adjusted OR per unit increase 1.53, 95% 
CI: 1.26, 1.86; P<0.001). The presence of cavitation on chest 
radiograph trended toward significance (P = 0.078), age of the 
contact (P = 0. 1 7) and SCG (p = 0.40) were not significandy 
associated with a HT household. We found no evidence of 
interaction between SCG and either cavitation (p = 0.50) or VACS 
(p = 0.49) but were underpowered for such comparisons and in 
fact had estimation problems. 

Discussion 

In this household contact study of smear-positive, culture 
confirmed pulmonary TB patients in Brazil, we systematically 
identified households with low and high transmission of MTB 
during a four-year period. Our study design and use of both TST 
and IGRA to measure prevalent and incident infection in exposed 
contacts showed heterogeneity of MTB transmission within 
households. We also found that cough severity in the index case 



is significandy associated with high transmission and that the 
presence of cavitation on chest radiograph is also associated, albeit 
marginally, with increased transmission. Our findings indicate that 
the infecting strain may result in distinct patterns of pathophys- 
iology of disease in index cases leading to differences in their 
infectiousness. 

The importance of cough frequency as a surrogate for TB 
infectiousness is well recognized [34,35,36,37]. Overnight cough 
frequency has been associated with the proportion of contacts with 
a positive TST [35]. However, the objective evaluation of cough, 
an inherendy subjective symptom, is challenging and we are 
unaware of previously published data relating cough strength or 
other measures of cough physiology with MTB transmission. Also, 
hosts may be more infectious during a specific time of day so there 
is a need for an objective and quantitative measure of infectious- 
ness [35,38]. Similarly, a hallmark of TB disease is the formation 
of lung cavitation, the result of liquefaction of tissue after the initial 
granulomatous response [39]. With progression of disease, 
cavitation becomes more common, occurring in 45-9 1 % of cases 
[40]; when present, cavities are associated with a higher bacillary 
load and increased transmission [35,39,41] when compared to 
patients with nodular lesions in the lung parenchyma [42,43,44]. 
However, measures of pulmonary TB disease severity such as 
cough strength, cavitation, and bacterial load in sputum are tighdy 
correlated and distinguishing their individual effects on MTB 
transmission can be challenging. We addressed this in our study 
design by using two separate instruments to measure cough 
(VACS and LCOJ, recording the presence of cavitations on chest 
radiograph and restricting enrollment to smear positive (—2+), 
culture-confirmed TB patients in an effort to maximize and 
homogenize already proven exposure across households. Whereas 
the VACS and LCQ^have been validated in other conditions with 
chronic cough (e.g. chronic obstructive pulmonary disease, 
asthma), to our knowledge this is the first study in which they 
have been applied to patients with pulmonary TB. 

There is accumulating evidence that strain diversity among 
MTB may result in clinically meaningful differences in disease 
outcome. Traditionally, MTB strain virulence has been defined 
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Figure 2. Cough severity (as measured by visual analog cough scale) in index TB case by household transmission category, single 
nucleotide polymorphism genetic cluster group (SCG), and cavitations on chest radiograph (N = 111). 

doi:1 0.1 371 /journal.pone.01 00984.g002 



based on phenotypic properties in varying conditions such as 
increased production of pro-inflammatory cytokines in the source 
or contact cases, growth characteristics in liquid cultures, 
macrophages, and THP1 cells, and growth kinetics in animal 
experiments [16,17,18]. Although some of the in-vitro data show 
conflicting results [45,46], current thinking is that virulent strains 
subvert the immune system by causing a reduced inflammatory 
response, which in turn results in impaired bacterial control by the 
host, more rapidly progressive disease, and enhanced transmission. 
Strain virulence also can be inferred from the capacity of MTB 
clones to spread in human populations, as measured by the 
identification of strain clusters using genotyping techniques 
[11,15,19,47]. However, it remains unclear whether there is a 
microbial basis to explain why some clinical isolates cause 



widespread disease and other closely related strains remain limited 
in spread [46,48,49]. Our study suggests that SCG 3b strains are 
associated with diminished transmission but that this may be due 
to weaker cough and decreased cavity formation, opening the 
possibility of strain-dependent variable lung parenchymal pathol- 
ogy and/or inflammation in the index TB patient. In fact, the 
hallmark of the local immune response in pulmonary TB is an 
unregulated activation of cells and expression of pro-inflammatory 
and anti-inflammatory cytokines and chemokines [50]. SCG 3b 
roughly corresponds to the H spoligotype; SCG5, the other SCG 
which was commonly observed in our index cases, roughly 
corresponds to the LAM spoligotype. Both LAM and H 
spoligo types are members of the Euro-American LSP/SNP based 
lineage [51]. This lineage is distinct from the East Asian 
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Table 2. Distribution of households according to study M. tuberculosis transmission categories by study time point. 





Study time point 




Final transmission category (2 nd TST) 




n (row %) 


Initial transmission category (1 st TST) 




Low 


Intermediate 


High 


n 


Low 


41 


23 (56) 


11 (27) 


7 (17) 


Intermediate 


30 


NA 


17 (57) 


13 (43) 


High 


53 


NA 


NA 


53 (100) 




Totals 


23 


28 


73 



NA = Not applicable. 

doi:1 0.1 371 /journal.pone.01 00984.t002 



tuberculosis lineages that includes Beijing type strains, which have 
been hypothesized to have increased virulence. Both SCG3b and 
SCG5 strains are expected to differ by a relatively small number of 
mutations because they are situated quite close to each other on 
SNP generated phylogenies, compared to more distantly related 
lineages such as the Beijing strains [26]. This suggests that only a 
small number of mutations may be responsible for the change in 
transmissibility observed in SC3b strains. These mutations may be 
discoverable in future whole genome sequencing studies. 

One important limitation has been that most studies examining 
the virulence potential of clustered MTB isolates have used banked 
isolates with limited knowledge about the epidemiologic and 
clinical conditions present at the time the strains were cultured, 
and thus fail to consider potential confounders related to 
characteristics of the index case (e.g. cough severity, bacterial 
load, cavitations), environment (e.g. ventilation, crowding) or 
contact (e.g. age, BCG and proximity of contact) [10]. Household 
contact studies such as this one are a classic, well-accepted design 



to study MTB transmission dynamics in a well-characterized 
group of individuals with a defined infectious challenge, while 
controlling for host, environmental and bacterial confounders. 
However, contacts that are already infected at the time of initial 
ascertainment may have been infected by the source case or, in an 
intermediate prevalence country such as Brazil, infected in the 
community [52,53]. As observed in our study, this translates into 
age-dependent differences in prevalence of infection in contacts 
(i.e. higher frequency of TB exposures in adults as a result of wider 
social networks). 

Our study has limitations. Whereas all patients were sputum 
smear- and culture-positive on solid media, we did not obtain 
quantitative culture data, which is known as a stronger surrogate 
of bacterial load and risk of infection. Between February 2008 and 
April 2009, we obtained a repeat TST in all contacts with TST < 
10 mm at study screening, including those living in IT households. 
As some IT households were excluded at the 1 st TST after April 
2009, some of which presumably would have been re-classified as 



Table 3. Household contacts with tuberculin skin test (TST) conversion from entry to week 8-12 by initial household transmission 
category and age of contact. 



TST conversion characteristic Initial (1" TST) household transmission category P* value 



No. with TST conversion/No. at risk (%) 




All 


Low 


Intermediate 


High 




Criterion I 1 


All 


50/271 (18.4) 


36/189 (19.0) 


13/72 (18.1) 


1/10 (9.1) 


0.67 


<5 years 


5/43 (11.6) 


4/32 (12.5) 


1/11 (9.1) 


0/0 (0.0) 


0.75 


6-15 years 


23/86 (26.7) 


17/63 (27.0) 


6/20 (30.0) 


0/3 (0.0) 


0.87 


5:16 years 


22/142 (15.5) 


15/94 (16.0) 


6/41 (14.6) 


1/7 (12.5) 


0.77 


Criterion 2 2 


All 


50/221 (22.6) 


36/1 59 (22.6) 


1 3/53 (24.5) 


1/9 (11.1) 


0.93 


<5 years 


5/42 (11.9) 


4/31 (12.9) 


1/11 (9.1) 


0/0 (0.0) 


0.72 


6-15 years 


23/76 (30.3) 


18/57 (31.6) 


5/16 (31.3) 


0/3 (0.0) 


0.87 


>16 years 


22/103 (21.4) 


14/61 (19.7) 


7/26 (26.9) 


1/6 (16.7) 


0.64 



n/N are number of converters divided by the number at risk as defined by the various criteria below. 

^ST conversion Criterion 1 (Brazilian guidelines}: 1 st TST<10 mm; 2 nd TST>10 mm; difference >10 mm. Does not include contacts with 1 st TST>10 mm (407, 56%), and 
those with missing 1 st (21, 3%) or 2 na TST (31, 4%) because they are not considered "at risk" by this criterion. 

2 TST conversion Criterion 2: 1 st TST<5 mm; 2 nd TST>10 mm; difference >6 mm. Does not include contacts with 1 st TST>5 mm (458, 65%), and those with missing 1 st 
TST {21, 3%) or 2 na TST (31, 4%) because they are not considered "at risk" by this criterion. 
*Score test from generalized estimating equation (GEE) estimation approach to logistic regression. 
Only contacts "at risk" of TST conversion are included for each criterion. 
doi:1 0.1 371 /journal.pone.01 00984.t003 
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HT, we cannot estimate the percentage of households that are HT 
within the study community. IT households are also underrepre- 
sented. However, because the initially HT group of households 
was large, the inclusion of additional households through 
conversion between the two TST applications into this group 
are unlikely to have altered the results of this study. Also, whereas 
our criteria to classify households were operationally driven, the 
definition of HT was robust in that 77% of HT households had 
100% of households with TST>10 mm. 

In conclusion, in this household contact study of culture 
confirmed, AFB 2+ smear-positive TB patients in a Brazilian 
urban setting we observed variable MTB transmission within 
households. Our results provide strong evidence of differential 
transmission among otherwise clinically similar patients with 
advanced TB disease. Data presented here also suggest that these 
differences in infectiousness may be modulated by yet unknown 
interactions between the host and the infecting strain, resulting in 
weaker cough and possibly, less formation of cavities. Further 
studies will be needed to mechanistically advance these novel 
findings. 

Supporting Information 

Table SI Additional characteristics of index tuberculo- 
sis cases, household contacts, and study dwellings at 
baseline by Mycobacterium tuberculosis transmission 
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